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We show how one can obtain a lower bound for the electrical, spin or heat conductivity of cor- 
related quantum systems described by Hamiltonians of the form H — Ho + gH\. Here Ho is an 
interacting Hamiltonian characterized by conservation laws which lead to an infinite conductivity 
for g = 0. The small perturbation gHi, however, renders the conductivity finite at finite temper- 
atures. For example, Hq could be a continuum field theory, where momentum is conserved, or an 
integrable one-dimensional model while Hi might describe the effects of weak disorder. In the limit 
g — > 0, we derive lower bounds for the relevant conductivities and show how they can be improved 
systematically using the memory matrix formalism. Furthermore, we discuss various applications 
and investigate under what conditions our lower bound may become exact. 
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I. INTRODUCTION 

Transport properties of complex materials are not only 
important for many applications but are also of funda- 
mental interest as their study can give insight into the 
nature of the relevant quasi particles and their interac- 
tions. 

Compared to thermodynamic quantities, the transport 
properties of interacting quantum systems are notori- 
ously difficult to calculate even in situations where in- 
teractions are weak. The reason is that conductivities 
of non-interacting systems are usually infinite even at fi- 
nite temperature, implying that even to lowest order in 
perturbation theory an infinite resummation of a per- 
turbative series is mandatory To lowest order this im- 
plies that one usually has to solve an integral equation, 
often written in terms of (quantum-) Boltzmann equa- 
tions or - within the Kubo formalism - in terms of vertex 
equations. The situation becomes even more difficult if 
the interactions are so strong that an expansion around 
a non-interacting system is not possible. Also numeri- 
cally, the calculation of zero-frequency conductivities of 
strongly interacting clean systems is a serious challenge 
and even for one-dimensional systems reliable calcula- 
tions are available for high temperatures onlyi2i£ii£i&. 

Variational estimates, e.g. for the ground state energy, 
are powerful theoretical techniques to obtain rigorous 
bounds on physical quantities. They can be used to 
guide approximation schemes to obtain simple analytic 
estimates and are sometimes the basis of sophisticated 
numerical methods like the density matrix renormaliza- 
tion groupi. 

Taking into account both the importance of transport 
quantities and the difficulties involved in their calculation 
it would be very useful to have general variational bounds 
for transport coefficients. 

A well known example where a bound for transport 
quantities has been derived is the variational solution 
of the Boltzmann equation, discussed extensively by 
Ziman 8 . The linearized Boltzmann equation in the pres- 



ence of a static electric field can be written in the form 
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where W^k' is the integral kernel describing the scatter- 
ing of quasiparticles and we have linearized the Boltz- 
mann equation around the Fermi (or Bose) distribution 

/£ = /°(ek) using / k = /° - £$ k . Therefore, the 

current is given by I = — e 2_, k v k ^-<I> k and the dc con- 
ductivity is determined from the inverse of the scattering 
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easy to see that this result can be obtained by maximiz- 
ing a functional^&i&li F[4>] with 
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where we used that W k , k / = reflecting the conser- 
vation of probability. The variational formula ((2J) is ac- 
tually closely related 8 - to the famous H-theorem of Boltz- 
mann which states that entropy always increases upon 
scattering. 

A lower bound for the conductivity can be obtained 
by varying $ only in a subspace of all possible func- 
tions. This allows for example to obtain analytically 
good estimates for conductivities without inverting an 
infinite dimensional matrix or, euqivalcntly, solving an 
integral equation, see Ziman's book for a large number 
of examples^. 

The applicability of Eq. © is restricted to situations 
where the Boltzmann equation is valid and bounds for 
the conductivity in more general setups are not known. 
However, for ballistic systems with infinite conductiv- 
ity it is possible to get a lower bound for the so-called 
Drude weight. Mazur— and later Suzuki^ considered 
situations where the presence of conservation laws pro- 
hibits the decay of certain correlation functions in the 
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FIG. 1: For g — a Drude peak shows up in the conductivity, 
resulting from exact conservation laws. For g the Drude 
peak broadens and the dc conductivity becomes finite. 

long time limit. In the context of transport theory their 
result can be applied to systems (see Appendix lAl for de- 
tails) where the finite-temperature conductivity a(ui,T) 
is infinite for to = and characterized by a finite Drude 
weight D{T) > with 

Re o{lo, T) = ttD{T)5(lu) + a Teg (u>, T). (3) 

Such a Drude weight can arise only in the presence of 
exact conservation laws Cj with [H, Cj] = 0. Suzuki 13 
showed that the Drude weight can be expressed as a sum 
over all Cj 

P ^ (C.J)2 p " {CjJ) 2 

where J is the current associated with a. For conve- 
nience a basis in the space of Ci has been chosen such 
that (CiCj) = for i ^ j. More useful than the equal- 
ity in Eq. (0| is often the inequality-^ which is obtained 
when the sum is restricted to a finite subset of conser- 
vation laws. Such a finite sum over simple expectation 
values can often be calculated rather easily using either 
analytical or numerical methods. The Mazur inequality 
has recently been used heavily^ i 14 ' 15 i 16 i 17 to discuss the 
transport properties of one-dimensional systems. 

Model systems, due to their simplicity, often exhibit 
symmetries not shared by real materials. For exam- 
ple, the heat conductivity of idealized one-dimensional 
Heisenberg chains is infinite at arbitrary temperature as 
the heat current is conserved. However, any additional 
coupling (next-nearest neighbor, inter-chain, disorder, 
phonon,...) renders the conductivity finite^^ ' 18 ! 19 ! 20 . 
If these perturbations are weak, the heat conductivity 
is, however, large as observed experimentally 21 ' 22 . For a 
more general example, consider an arbitrary translation- 
ally invariant continuum field theory. Here momentum 
is conserved which usually implies that the conductivity 
is infinite for this model. In real materials momentum 
decays by Umklapp scattering or disorder rendering the 
conductivity finite. It is obviously desirable to have a 
reliable method to calculate transport in such situations. 
In this work we consider systems with the Hamiltonian 



where for g = the relevant heat-, charge- or spin- con- 
ductivity is infinite and characterized by a finite Drude 
weight given by Eq. Q. As discussed above, H might be 
an integrable one-dimensional model, a continuum field 
theory, or just a non-interacting system. The term gH\ 
describes a (weak) perturbation which renders the con- 
ductivity finite, e.g. due to Umklapp scattering or dis- 
order, see Fig. [TJ Our goal is to find a variational lower 
bound for conductivities in the spirit of Eq. (J5J) for this 
very general situation, without any requirement on the 
existence of quasi particles. For technical reasons (see 
below) we restrict our analysis to situations where H is 
time reversal invariant. 

In the following, we first describe the general setup and 
introduce the memory matrix formalism, which allows 
us to formulate an inequality for transport coefficients 
for weakly perturbed systems. We will argue that the 
inequality is valid under the conditions which we specify. 
Finally, we investigate under which conditions the lower 
bounds become exact and briefly discuss applications of 
our formula. 



II. SETUP 

Consider the local density p{x) of an arbitrary phys- 
ical quantity which is locally conserved, thus obeying a 
continuity equation 

d t p + Vj = 0. 

Transport of that quantity is described by the dc con- 
ductivity a which is the response of the current to some 
external field E coupling to the current, 

(J) = aE, 

where J = J j(x) is the total current and (J) its expec- 
tation value. Note that J can be an electrical-, spin-, or 
heat current and E the corresponding conjugate field de- 
pending on the context. The dynamic conductivity <j{z) 
is given by Kubo's formula, see Eq. (|A1|) . We are inter- 
ested in the dc conductivity a = lim^^o f {z — lj + iO). 

Starting from the Hamiltonian (J5J) we consider a sys- 
tem where Hq posesses a set of exact conservation laws 
{Ci} of which at least one correlates with the current, 
(JCi) 7^ 0. Without loss of generality we assume 
(CiCj) = for i f j. For g = the Drude weight D 
defined by Eq. §5§ is given by Eq. (H]). We can split up 
the current under consideration into a part which is par- 
allel to the Ci and one that is orthogonal, 

J = J\\ + J±, 

with J|| = 2^ \q2\ Ci, which results in a separation of 
the conductivity, 



H = H +gH 1 , 



(5) 



a(z) = a\\(z) +<r±(z). 



(6) 
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Since the conductivity a(z) is given by a current-current 
correlation function and the current Jh (J±) is diago- 
nal (off-diagonal) in energy, cross-correlation functions 
((Jm; J_l)) vanish in Eq. ©. 

According to Eq. ((4]), the Drude peak of the unper- 
turbed system, g = 0, arises solely from J||: 



Re CT|| (w) = ttD5(lo), 



(7) 



while crj^z) appears in Eq. ([3]) as the regular part, 
Re cr±(u>) — a Icg (uj). 

In this work we will focus on a\\{uj), since the small 
perturbation is not going to affect <7j_(w) much (which is 
assumed to be free of singularities here, see section IIV[) 
while cry (to = 0) diverges for g — > 0, see Fig. [T] As we 
are interested in the small g asymptotics only, we may 
neglect the contribution cr±(0) to the dc conductivity. 
Hence we set J = Jh and <j{oj) = <J\\ (lo) in the following. 



III. MEMORY MATRIX FORMALISM 

We have seen that certain conservation laws of Ho play 
a crucial role in determining the conductivity of both the 
unperturbed and the perturbed system. In the presence 
of a small perturbation gH\, these modes are not con- 
served anymore but at least some of them decay slowly. 
Typically, the conductivity of the perturbed system will 
be determined by the dynamics of these slow modes. To 
separate the dynamics of the slow modes from the rest, it 
is convenient to use a hydrodynamic approach based on 
the projection of the dynamics onto these slow modes. In 
this section we will therefore review the so called memory 
matrix formalism 2 ^, introduced by Mori and Zwanzig 24-25 
for this purpose. In the next section we will show that 
this approach can be used to obtain a lower bound for 
the dc conductivity for small g. 

We start by defining a scalar product in the space of 
quantum mechanical operators, 



(A\B) 



dX{A^B(i\)) - I3{A^){B) 



(8) 



As the next step we choose a - for the moment - arbi- 
trary set of operators {Ci}. In most applications, the Ci 
are the relevant slow modes of the system. For notational 
convenience, we assume that the {Ci} are orthonormal- 
ized, 



(Ci\Cj) = 5. K 



(9) 



In terms of these we may define the projector P onto (and 
Q away from) the space spanned by these 'slow' modes 

P = J2\Ci)(Ci\ = l-Q. 



We assume that C\ is the current we are interested in, 
| J) = \C 1 ). 



The time evolution is given by the Liouville- 
(super)operator 

L = [H, .] = L + gh x 

with (LA\B) = (A\LB) = (A\L\B), and the time evo- 
lution of an operator may be expressed as |^4(i)) = 



Ae 



-iHt\ 



e \A). With these notions, one obtains 



the following simple, yet formal expression for the con- 
ductivity: 



J 



L 



J 



Using a number of simple manipulations, one can 
sho w 23 ' 24 i 25 that the conductivity can be expressed as 
the (1, l)-component of the inverse of a matrix 



cr(w) = (M(uj) + iK — ioj) n , 



where 



M«(w)= C, 



' lu — LQ 



C, 



is the so-called memory matrix and 



(10) 



(11) 



(12) 



a frequency independent matrix. The formal expression 
(JTOJ) for the conductivity is exact, and completely gen- 
eral, i.e. valid for an arbitrary choice of the modes Ci 
(they do not even have to be 'slow'). Only C\ = J is re- 
quired. However, due to the projection operators Q, the 
memory matrix (jlip is in general difficult to evaluate. It 
is when one uses approximations to M that the choice of 
the projectors becomes crucial (see below). 

Obviously, the dc conductivity is given by the (1, 1)- 
component of 



(M(0) + K)-\ 



(13) 



More generally, the (m, n)-component of Eq. (|13p de- 
scribes the response of the 'current' C m to an external 
field coupling solely to C n . We note, that since a matrix 
of transport coefficients has to be positive (semi)definitc, 
this also holds for the matrix M(0) + K. 

To avoid technical complications associated with the 
presence of K we restrict our analysis in the following to 
time reversal invariant systems and choose the Ci such, 
that they have either signature +1 or —1 under time 
reversal^ O. In the dc limit, u> = 0, components of 
Eq. (fT"3"|) connecting modes of different signatures vanish. 
Thus, M(0) + K is block-diagonal with respect to the 
time reversal signature, and consequently we can restrict 
our analysis to the subspace of slow modes with the same 
signature as C\. However, if C m and C n have the same 
signature, then (C m \C n ) = 0, and thus K vanishes on 
this restricted space. The dc conductivity therefore takes 
the form 



(M(O)- 1 ) 



li • 



(14) 
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IV. CENTRAL CONJECTURE 

To obtain a controlled approximation to the memory 
matrix in the limit of small g, it is important to identify 
the relevant slow modes of the system. For the C, we 
choose quantities which are conserved by Hq, [Ho, C»] = 
0, such that C7j = ig[Hi, Cj] is linear in the small coupling 
g. As argued below, we require that the singularities 
of correlation functions of the unperturbed system are 
exclusively due to exact conservation laws Ci, i.e. that 
the Drude peak appearing in Eq. |3]) is the only singular 
contribution. Furthermore, we choose J = Ju = C\ and 
consider only C; with the same time reversal signature 
as J, as discussed in the previous section. 

To formulate our central conjecture we introduce the 
following notions. We define M n (ui) as the (exact) n x n 
memory matrix obtained by setting up the memory ma- 
trix formalism for the first n slow modes {Cj, i = 1, ..,n}. 
Note that the definitions of the relevant projectors P and 
Q also depend on this choice, and that for any choice of 
n one gets a — (M~ 1 )u. We now introduce the approxi- 
mate memory matrix M n motivated by the following ar- 
guments: Ci is already linear in g, therefore in Eq. (jlip 
we approximate L by Lq and replace (.|.) by (.|.)o as 
we evaluate the scalar product with respect to Hq. As 
Lo\Ci) = and (Cj\Ci) = due to time reversal symme- 
try, one has LqQ = 1 and Q\Ci) = \Ci) and therefore the 
projector Q does not contribute within this approxima- 
tion. We thus define the n x n matrix M n by 



M„ 



lim C 



U) - Lq 



c, 



(15) 



Note that M n is a submatrix of M m for m > n and 
therefore the approximate expression for the conductiv- 
ity a (M~ 1 )u does depend on n while (M~ 1 )u is 
independent of n. A much simpler, alternative deriva- 
tion for Mi is given in Appendix [B] where the validity of 
this formula is also discussed. 

The central conjecture of our paper is, that for small g 
(M~ 1 )n gives a lower bound to the dc conductivity, or, 
more precisely, 



\l/g* 



(M- 1 ) 11 >--->(Mn 1 )ll> 



> Mf 1 . 



(16) 



Here a 



1/9 2 



(1/g 2 ) lim g ^ g 2 c denotes the leading 
term oc l/<? 2 in the small-g expansion of a. Note that 
M n oc g 2 by construction. M^ is the approximate mem- 
ory matrix where al$^ conservation laws have been in- 
cluded. In some special situations, discussed in Ref. 0, 
one has a ~ 1/g 4 and therefore &\i/ g 2 — oo. 

A special case of the inequality above is Eq. (|B4|) in 
appedix [Bl as the scattering rate T/x may be expressed 
as f/x 2 = Mi. 

Two steps are necessary to prove Eq. (TT^l) . The simple 
part is actually the inequalities in Eq. (flu]) . They are 
a consequence of the fact that the matrices M n are all 



positive definite and that M n is a submatrix of M m for 
m > n. More difficult to prove is that the first equality 
in (TIB|) holds. To show this we will need an additional 
assumption, namely, that the regular part of all correla- 
tion functions (to be defined below) remains finite in the 
limit g — > 0, u> — > 0. In this case, the perturbative ex- 
pansion around M^ in powers of g is free of singularities 
at finite temperature (which is not the case for M n<oc ). 
This in turn implies that limg_>o M^/g 2 = M^lg 2 and 
therefore cr| 1/g2 = {M^)\\. 

Next, we present the two parts of the proof. 



A. Inequalities 

We start by investigating the (l,l)-component of the 
inverse of the positive definite symmetric matrix M^ . It 
is convenient to write the inverse as 



(M- l )u = max % Gl) 

V If M ooV 



(17) 



where e\ is the first unit vector. The same method is 
used to derive Eq. ^ in the context of the Boltzmann 
equation. The maximum is obtained for ip = M^e\. 
By restricting the variational space in (jTTJ) to the first n 
components of ip we reproduce the submatrix M n of M^ 
and obtain 



> 



max ^ 11 = (M~ )n 

~FTx = i M n<m)n 

<p=T,i K Vie* <P Moop 



By choosing different values for m and n < m, this proves 
all inequalities appearing in (p~6|) . 



B. Expansion of the memory matrix 

We proceed by expanding the exact memory matrix 
M n , where P n = 1 — Q n is a projector on the first n 
conservation laws, in powers of g. Using that LQ n = 
Lq + gLiQ ni we obtain the geometric series 



fc=0 



Ci Qn— ~j— ( LiQ„— — ) Cj 



Lr 



Note that this is not a full expansion in g, as the scalar 
product ([5]) is defined with respect to the full Hamilto- 
nian H = Hq + gH\. We will turn to the discussion of 
the remaining ^-dependence later. 
In general, one can expand 



L\ — ^ ' A TOn |j4 m ) (A n | 
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in terms of some basis A m in the space of operators. 
Therefore Eq. (TT5)) can be written as a sum over products 
of terms with the general structure 



1 



U — Lq 



D 



(19) 



In the following we would like to argue that such an ex- 
pansion is regular for n — oo if all conservation laws 
have been included in the definition of Q. As argued in 
Appendix[Bl we have to investigate whether the series co- 
efficients in Eq. (fT5|) diverge for to — > 0. The basis of our 
argument is the following: as Qoo projects the dynam- 
ics to the space perpendicular to all of the conservation 
laws, the associated singularities are absent in Eq. (TlT)|) 
and therefore the expansion of is regular. 

To show this more formally, we split up B = B\\ + B± 
in (|19[) into a component parallel and one perpendicular 
to the space of all conserved quantities, \Bn) — P^B). 
With this notation, the action of Lq becomes more trans- 
parent: 



1 



u - L n 



\B) 



1 



1 



\B i 



(20) 



As we assume that all divergencies can be traced back 
to the conservation laws, we take the second term to be 
regular. It is only the first term which leads in Eq. (fT9")) 
to a divergence for to — > 0, provided that (A|<5„|B|i) is fi- 
nite. If we consider the perturbative expansion of M n<oc , 
where P n = 1 — Q n projects only to a subset of con- 
served quantities, then finite contributions of the form 
(yl|<5„|i?|| ) exist and the perturbative series in g will be 
singular (see also Appendix [BJ. Considering Moo, how- 
ever, Qoo projects out all conservation laws and therefore 
by construction Qoc\B\\) = QooPoo\B) — 0. Thus the 
first term in (|20f does not contribute in (fT9j) for n — oo 
and the expansion (TTg)) of is therefore regular. 

The only remaining part of our argument is to show 
that in the limit g — > one can safely replace (.|.) by 
(.|.)o- Here it is useful to realize that {A\B) can be in- 
terpreted as a (generalized) static susceptibility. In the 
absence of a phase transition and at finite temperatures, 
susceptibilities are smooth, non-singular functions of the 
coupling constants and therefore we do not expect any 
further singularities from this step. If we define a phase 
transition by a singularity in some generalized suscepti- 
bility, then the statement that susceptibilities are regular 
in the absence of phase transitions even becomes a mere 
tautology. 

Combining all arguments, the expansion (TT5)) of 
— > 0) is regular, and using {Gi\Q oo — (C{\ [see 
discussion before Eq. (fT5]) ] its leading term, k = is 
given by . We therefore have shown the missing first 
equality of our central conjecture (|16p . 



V. DISCUSSION 

In this paper we have established that in the limit of 
small perturbations, H — Hq + gH\, lower bounds to dc 
conductivities may be calculated for situations where the 
conductivity is infinite for g = 0. In the opposite case, 
when the conductivity is finite for g — 0, one can use 
naive perturbation theory to calculate small corrections 
to a without further complications. 

The relevant lower bounds are directly obtained from 
the memory matrix formalism. Typicall y 26 i 27 i 28 one has 
to evaluate a small number of correlation functions and 
to invert small matrices. The quality of the lower bounds 
depends decisively on whether one has been able to iden- 
tify the 'slowest' modes in the system. 

There are many possible applications for the results 
presented in this paper. The mostly considered situ- 
ation is the case where Hq describes a non-interacting 
system^. For situations where the Boltzmann equation 
can be applied, it has been pointed out a long time ago by 
Belitz 2 - that there is a one-to-one relation of the memory 
matrix calculation to a certain variational Ansatz to the 
Boltzmann equation, see Eq. @. In this paper we were 
able to generalize this result to cases where a Boltzmann 
description is not possible. For example, if Hq is the 
Hamiltonian of a Luttinger liquid, i.e. a non-interacting 
bosonic system, then typical perturbations are of the 
form cos <f> for which a simple transport theory in the 
spirit of a Boltzmann or vertex equation does not exist 
to our knowledge. 

Another class of applications are systems where Hq 
describes an interacting system, e.g. an integrable one- 
dimensional model^ or some non-trivial quantum-field 
theory^. In these cases it can become difficult to cal- 
culate the memory matrix and one has to resort to use 
either numerical^ or field-theoretical methods^ to obtain 
the relevant correlation functions. 

An important special case are situations where Hq 
is characterized by a single conserved current with the 
proper symmetries, i.e. with overlap to the (heat-, spin- 
or charge-) current J. For example, in a non-trivial con- 
tinuum field theory Hq, interactions lead to the decay of 
all modes with exception of the momentum P. In this 
case the momentum relaxation and therefore the con- 
ductivity at finite T is determined by small perturba- 
tions gH\ like disorder or Umklapp scattering which are 
present in almost any realistic system. As = M\ in 
this case, our results suggest that for small g the conduc- 
tivity is exactly determined by the momentum relaxation 
rate M PP = linwo i(P\(uj - Lq)- 1 ^), 



Xpj 
M P P 



for g^0. 



(21) 



Here we used that Jy = P(P\J)/(P\P) with X pj = (P\J) 
and we have restored all factors which arise if the nor- 
malization condition ([5]) is not used. In Appendix [Cl we 
check numerically that this statement is valid for a real- 
istic example within the Boltzmann equation approach. 
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A number of assumptions entered our arguments. The 
strongest one is the restriction that all relevant singu- 
larities arise from exact conservation laws of Ho. We 
assumed that the regular parts of correlation functions 
are finite for lo = 0. There are two distinct scenarios 
in which this assumption does not hold. First, in the 
limit T — > 0, often scattering rates vanish which can lead 
to diverencies of the nominally regular parts of correla- 
tion functions. Furthermore, at T = even infmitesi- 
mally small perturbations can induce phase transitions 
- again a situation where our arguments fail. Therefore 
our results are not applicable at T = 0. Second, finite 
temperature transport may be plagued by additional di- 
vergencies for uj — ► not captured by the Drude weight. 
In some special models, for instance, transport is singu- 
lar even in the absence of exactly conserved quantities 
(e.g. non-interacting phonons in a disordered crystal 8 ). 
In all cases known to us, these divergencies can be traced 
back to the presence of some slow modes in the system 
(e.g. phonons with very low momentum). While we have 
not kept track of such divergencies in our arguments, we 
nevertheless believe that they do not invalidate our main 
inequality (|16J) as further slow modes not captured by ex- 
act conservation laws will only increase the conductivity. 
It is, however, likely that the equality pip is not valid 
for such situations. In Appendix [C] we analyze in some 
detail within the Boltzmann equation formalism under 
which conditions (f2Tj) holds. As an aside, we note that 
the singular heat transport of non-interacting disordered 
phonons, mentioned above, is well described within our 
formalism if we model the clean system by Hq and the 
disorder by H\, see the extensive discussion by Ziman 8 
within the variational approach which can be dir ectly 
translated to the memory matrix language, see Ref. [291 ] . 

It would be interesting to generalize our results to cases 
where time reversal symmetry is broken, e.g. by an exter- 
nal magnetic field. As time reversal invariance entered 
nontrivially in our arguments, this seems not to be sim- 
ple. We nevertheless do not see any physical reason why 
the inequality should not be valid in this case, too. One 
example where no problems arise are spin chains in a 
uniform magnetic field— where one can map the field to 
a chemical potential using a Jordan- Wigner transforma- 
tion. Then one can directly apply our results to the time 
reversal invariant system of Jordan- Wigner fermions. 



APPENDIX A: DRUDE WEIGHT AND MAZUR 
INEQUALITY 

In this appendix we clarify the connection between the 
Drude weight and the Mazur inequality, mentioned in the 
introduction. 

The Drude weight D is the singular part of the con- 
ductivity at zero frequency, Re cr(u>) = ttD6(lj) +a ies (u)). 
It can be calculated from the relation 

D = lim lu Im ct(lu). 

It has been introduced by Kohn 32 as a measure of ballistic 
transport, indicated by D > 0. 

Using Kubo formulas, conductivities can be expressed 
in terms of the dynamic current susceptibilities^ 3 - Tl(z) 
using 

a(z) = -l(n T -n(z)), (Al) 
where II(z) is the current response function 

R(z) = -jf dte^([J(t),J(0)}) (A2) 

rF = /f£n>). (A3) 

J 7T LU 

and II T is a current susceptibility. The conductivity may 
be calculated by setting a(ui) — a{z — lu + iO). Relation 
(|A3p is a well known sum rule and for all regular corre- 
lation functions one has n T = 11(0). In the presence of 
a singular contribution to o~{lu), one easily identifies the 
Drude weight with the expression II T — 11(0). For this dif- 
ference Mazu r 12 i 13 derived a lower bound. Furthermore, 
Suzukii 3 . has shown, that II T - 11(0) may be expressed 
as a sum over all constants of the motion Cj present in 
the system 3 -^, 

^ = n T -n ( o)^f;^. (A4) 

n=0 ^ i' 

Thus, the Drude weight is intimately connected to the 
presence of conservation laws: only components of the 
current perpendicular to all conservation laws decay and 
any conservation law with a component parallel to the 
current (i.e. with a finite cross-correlation (CjJ)) leads 
to a finite Drude weight and thus ballistic transport. The 
relation between the Drude weight and Mazur 's inequal- 
ity has been first pointed out by ZotosJ^. 
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APPENDIX B: PERTURBATION THEORY FOR 

l/o 

Let us give an example of a naive perturbative deriva- 
tion (see also Ref. @) to gain some insight about what 
problems can turn up in a perturbative derivation as the 
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one presented in this work. According to our assump- 
tions, the conductivity is diverging for g — > and there- 
fore it is useful to consider the scattering rate T{w)/\ 
(with the current susceptibility x) defined by 



ct(w) 



X 



T(uj)/x-iu' 



(Bl) 



If J is conserved for g = (i.e. for J = Jy, see above), 
the scattering rate vanishes, T(uj) = 0, for g = 0, which 
results in a finite Drude weight. A perturbation around 
this singular point results in a finite T(uj). In the limit 
g — > we can expand (|Bip for any finite frequency a; in 
r to obtain 



u 2 Re cr(w) = Re T(w) + 0{T 2 /uj). 



(B2) 



We can read this as an equation for the leading order 
contribution to r(cj), which now is expressed through 
the Kubo formula for the conductivity. By partially in- 
tegrating twice in time we can write r(a;) = T(uj) + 0(g 3 ) 
with 



Ref(w) = Re - — 
z V 



dte tzt ([j(t), j(0)]) 

z— oj-\-iQ 

(B3) 

where J = i[H, J) = ig[H±, J] is linear in g and therefore 
the expectation value (...)o can be evaluated with respect 
to Hq (which may describe an interacting system). Thus 
we have expressed the scattering rate via a simple corre- 
lation function of the time derivative of the current. 

To determine the dc conductivity one is interested in 
the limit u> — > and it is tempting to set u> = in 
Eq. (|B3p . We have, however, derived Eq. (|B3p in the 
limit g — > at finite u> and not in the limit uj — > at 
finite g. The series Eq. (|B2[) is well defined for finite 
lu =/= only and in the limit to —> the series shows 
singularities to arbitrarily high orders in 1/lo. 

At first sight this makes Eq. (|B3[) useless for calculating 
the dc conductivity. One of the main results of this paper 
is that, nevertheless, T(uj = 0) can be used to obtain a 
lower bound to the dc conductivity 



a(w = 0) > W 



X 



r(o) 



for g — > 0. 



(B4) 



APPENDIX C: SINGLE SLOW MODE 

In this appendix we check whether in the presence of 
a single conservation law with finite cross correlations 
with the current the inequality (|16p can be replaced by 
the equality (|2ip . This requires us to compare the true 
conductivity, which in general is hard to determine, to 
the result given by M\. Thus we restrict ourselves to 
the discussion of models for which a Boltzmann equation 
can be formulated and the expression for the conductivity 
can be calculated at least numerically. In the following 
we first show numerically that the equality ([2"T]) holds for 



a realistic model. In a second step we discuss the precise 
regularity requirement of the scattering matrix such that 
Eq. (SI} holds. 

To simplify numerics, we consider a simple one- 
dimensional Boltzmann equation of interacting and 
weakly disordered Fermions. Clearly, the Boltzmann ap- 
proach breaks down close to the Fermi surface due to 
singularities associated with the formation of a Luttinger 
liquid, but in the present context we are not interested 
in this physics as we only want to investigate properties 
of the Boltzmann equation. To avoid the restrictions as- 
sociated with momentum and energy conservation in one 
dimension we consider a dispersion with two minima and 
four Fermi points, 



k 2 k 4 1 



a- 



2 4 10 
The Boltzmann equation reads 

df° k 



(CI) 



dek 



tf=£^[/*/v(i-/ 9 )(i-/ 9 o 



k'qq' 



-/,/,'(! -/*)(! -A-')] 
■ g 2 S(e k - e k ,)[f k (l - f k ,) - fa (1 - fa 



(C2) 



where the inelastic scattering term Sf^, = 5(e k + e k ' — 
e q — e q /)5(k + k' — q — q') conserves both energy and mo- 
mentum. In the last line we have linearized the right 
hand side using the definitions of the introductory chap- 
ter. The velocity v k is given by Vk = -^f-k- The scat- 
tering matrix splits up into an interaction component 
and a disorder component, W k k> = W^ k , + g 2 W k 1 k ,. As 
we do not consider Umklapp scattering, W kk i conserves 
momentum, J^k' W kk /k' = 0, and one expects that mo- 
mentum relaxation will determine the conductivity for 
small g. 

For the numerical calculation we discretize momentum 
in the interval [— ir/2, 7r/2], k n — nSk — nn/N with inte- 
ger n. (At the boundaries the energy is already too high 
to play any role in transport.) The delta function aris- 
ing from energy conservation is replaced by a gaussian of 
width S. The proper thermodynamic limit can for exam- 
ple be obtained by choosing S — 0.3/yN. The numerics 
shows small finite size effects. 

In Fig. [5] we compare the numerical solution of the 
Boltzmann equation to the single mode memory matrix 
calculation or, equivalently 29 , to the variational bound 
obtained by setting = k in Eq. @ 



k K de k 



Ek.k'kW kk ,k' g 2 Tk.k-kWl k ,k'- 



(C3) 



As can be seen from the inset, in the limit of small g 
one obtains the exact value for the conductivity, which is 
what we intended to demonstrate. 
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under the condition that 



FIG. 2: Comparison of the result of a single mode memory 
matrix calculation (solid line), Eq. (|C3f) . to the full numerical 
solution of the Boltzmann equation (dotted line) for T = 0.05 
and N = 500. The memory matrix is always a lower bound to 
the Boltzmann result and converges towards it as the disorder 
strength g is reduced, as shown in the inset (ratio of the single 
mode approximation to the Boltzmann result). 



Next we turn to an analysis of regularity conditions 
which have to be met in general by the scattering matrix 
Wkk' such that convergence is guaranteed in the limit 
g — ► 0. According to the assumptions of this appendix, 
for g — the variational form of the Boltzmann equa- 
tion ([2]) has a unique solution $fc (up to a multiplica- 
tive constant), with F($ k ) = oo, J2k< W kk>®k> = and 
£ fc i>k**4f%fe*>0. 

In the presence of a finite, but small g we write the solu- 
tion of the Boltzmann equation as $ = $+$ , where $ ± 
has no component parallel to $ (i.e. J^k ^k^kdf / dek = 
). On the basis of the two inequalities 



F[<f>] < 



$ i 14/$ 



(C4) 

> QgZW 1 ® (C5) 



one concludes that Eq. (|2ip is valid, i.e. that 

F[$] 



lim 



1 



,. * 4f i df° 

hm > w fc $fe— = > vA-f- 

9 ^°k dek k dek 



or, equivalcntly, 



o 



,. * d f 

hm ) « t $y — 

9^° ^ de k 



(C6) 



We therefore have to check whether Q± becomes small 
in the limit of small g. 

Expanding the saddlepoint equation for ([2]) we obtain 



As by definition has no component parallel to <E>, we 
can insert the projector Q which projects out the con- 
servation law in front of $^ on the left hand side. We 
therefore conclude that if the inverse of W°Q exists, then 
<f>j^ is of order g 2 , Eq. (|C6[) is valid and therefore also 
Eq. (|21|) . In our numerical examples these conditions are 
all met. 



Under what conditions can one expect that Eq. ()C6|) is 
not valid? Within the assumptions of this appendix we 
have excluded the presence of other zero modes of W° 
(i.e. conservation laws) with finite overlap with the cur- 
rent. But it may happen that W° has many eigenvalues 
which are arbitrarily small such that the sum in Eq. (|C6[) 
diverges. In such a situation the presence of slow modes 
which cannot be identified with conservation laws of the 
unperturbed system invalidates Eq. (|2ip . 
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